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UNIT VI: THE LINEAR SYSTEMS 



We study linear systems of n first-order differential equations. They are related to first-order 
matrix differential equations. When the corresponding matrix is constant, then the eigenvalues 
and the eigenfunctions of the matrix provide a useful framework to construct the general solution. 
The fundamental matrix is constructed as the exponential matrix. 



LECTURE 25. LINEAR SYSTEMS 

A linear system of n differential equations in n unknowns is 

y'l = aii{t)yi H h ai„(t)y„ + 

2/2 = «2l(t)yi H h a2n{t)yn + f2{t), 

(25.1) 

y'n = anl{t)yi H h ann{t)yn + fn{t). 

Here and elsewhere, ' = d/dt. In the matrix notation, 

(25.2) y' = A{t)y + f{t), 

where y = . . . , yn{t))'^ and f{t) = {fi{t), . . . , fn{t))'^ are vector-valued functions defined 

on an interval t G / and taking values in M", and A{t) = j=i is an n x n matrix-valued 

function defined on /. 

A matrix-valued function A{t) = {aij{t)) is said to be continuous, bounded, or differentiable if 
each element aij of A{t) is continuous, bounded and differentiable, respectively. Differentiation 
and integration are element- wise: 

Let us define an operator 

(25.3) Ly = y'-Ay. 

With this notation, (25.2) is written as Ly = f{t). The domain of the operator L is the space of 
n-dimensional vector-valued functions differentiable on /. 

Exercise. Show that L is linear. That is, 

L{ciyi{t) + C2y2{t)) = ciLyi{t) + C2Ly2{t). 

Since L is linear, the fundamental principle of superposition and the principle of the comple- 
mentary function and other results for linear operators apply. 

Existence and Uniqueness result. If A{t) and f{t) are continuous and bounded on an interval I 
containing to, then for each ijQ G M" the initial value problem 

y' = A{t)y + f{t), y{to) = yo 
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has a unique solution in /. 

Working assumption. Throughout the note and the following notes, A{t) and f{t) are assumed 
to be continuous and bounded on an interval tel. 

Linear independence. The vecotrs yi , . . . , y„ in M" are said to be linearly independent if 

ciyi + C2y2 H h CnVn = implies ci = C2 = • • • = c„ = 0. 

Let Y = (yi, . . . , yn) be an n X n matrix whose j-th column is yj. Then, yi, . . . , y„ are linearly 
independent if and only if the determinant |y| / 0. In this case, moreover, these functios form a 
basis for the linear space M". 

Similarly, the vector-valued functions yi{t), . . . ,yn{t) are said to be linearly independent on the 
interval / if 

ciyi(i) + C2?/2(i) H hcnyn(t) = on / implies ci = C2 = • • • = c„ = 0. 

This definition is more restrictive than that for vectors. But, if the functions are solutions of linear 
systems of differential equations, then their linear independence can be characterized in terms of 
the determinant of a matrix-valued function, analogous to that for vectors. 

A basis of solutions of 



where A{t) is a continuous n x n matrix-valued function on /, is a collection of n solutions 
yi(t), . . . , yn{t) so that every solution can be expressed uniquely as their linear combination 



for constants ci, C2, . . . , c„. 

Theorem 25.1. Let yi{t), . . . , y„(t) be n solutions of (25.4) on I. Let Y{t) = {yi{t), . . . , yn{t)) be the 
n X n matrix-valued function, whose j-th column is yj{t). 

The functions form a basis of solutions if and only z/|y(to)| ^ at some point to G /. If\Y{tQ)\ = at 
a point to G I, then \Y{t)\ = Ofor all t £ I. 

If |l^(to)| = then the column vectors in l"(to) are linearly dependent. Thus, we can find 
constants cj not all zero such that 



satisfies y(to) = 0. Moreover, y{t) is a solution of (25.4). Then, by uniqueness, y{t) = for all tel. 
That is, the functions are linearly dependent on /. 

Notations. \Y{t)\ = det(yi(t) • • • y„(t)) is called the Wronskian of yi{t), . . . , yn{t)- 

If the columns of Y{t) are linearly independent solutions of y ' = A{t)y then Y{t) is called the 
fundamental matrix for y ' = A{t)y. The fundamental matrix yields the general solution, and hence 
the complete solution, of y ' = A{t)y. This proves the assertion. 

Plane autonomous systems. When n = 2, then the system 

y'l = aii(t)yi + ai2{t)y2, 
y2 = a2i{t)yi + a22{t)y2 

is called a plane system. In the matrix notation, we write it as y ' = A{t)y, where y = (yi, y2) and 



(25.4) 



y = Ait)y, 



y{t) = ciyi(t) H h c„yn(t) 



y{t) = ciyi(t) H h c„yn(t) 
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When Uij are constants, then let 4> = yi. We eHminate 2/2 from the system, and 

{(j) - aii4>)' = 012021'/' + 022(0' - aii4>)- 
Thus, we obtain the companion or secular equation 
(25.6) 0" - (tr A)<j)' + (det A)0 = 

of (25.5). 

Proposition 25.2. If {yi{t) , y2{t))'^ is a solution of (25.5), then both yi{t) and y2{t) are solutions of the 
secular equation (25.6). 

The proof is easy and it is left as an exercise. 

Conversely, the secular equation (25.6) can be used to solve the system (25.5), by finding the 
roots of the corresponding characteristic polynomial 

- (trA)A + det^, 

and subsequently (25.5) can be solved. In this note, we caU the quadratic polynomial as the char- 
acteristic polynomial of A. 

Example 25.3. We consider the linear system 

y'l = 12yi +5y2, ^ _ /12 5 

2/2 = -%! + 2/2, ^ 

Its secular equations is (j)" — 130' + 420 = 0, and the solutions are generated by e^* and e''*. 

If yi{t) = cie^*, then 5y2 = y'l — 12yi = — 6cie^*. If yi{t) = C2e'^* then y2{t) = — C2e''*. Thus, the 
general solutions of the system is 



Exercise. Find the general solution of 

1. A= (^^^ . (Answer, cie* Q + cae* (^_^*^ ^ ). 

2- ^ = (-4 1) • (^^S"^"^- ^1^* (TsinL) + {2Zsl) 

Review of linear algebra. We end this note by listing notations and results of linear algebra, 
useful in the study of linear system of differential equations. 

Let A = {aij) be an n x n matrix. If the rows and the columns of A are interchanged then the 
resulting matrix is called the transpose of A, denoted by A^. It is immediate that 

{A^f = A, (^^)-i = {A-y, {ABf = B^A'. 

A matrix is said to be symmetric if ^4 = A^. 

We use \ A\ for the determinant of A. The {i, j)-minor of A, denoted by Mij, is the determinant 
of the (n — 1) X (n — 1) matrix obtained when the ith row and the jth column of A are deleted. 
Then, by definition, 

\A\ = auMu - ai2Mi2 + • • • + (-l)"+^ai„Mi„. 
The cofactor of A, denoted by cof A,isanxn matrix {cij)fj^^, where 

= i-iy+'Mij. 
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Then, 

1^1 = OllCll + ai2Ci2 H h ainCin- 

Finally, the adjugate of A, denoted by adj A, is the transpose of the cofactor matrix 

adj A = (cof A)"^, (adj A)ij = cji. 
Then, follows the Laplaces expansion formula for determinant 

n n 

\A\ = ^aij{cofA)ij = ^aij{cofA)ij 
j=i 1=1 

for any 1 ^ i, j ^ n. That is a determinant can be computed on any row or on any column. 
Example 25.4. The adjugate of the 2x2 matrix A = !/ ^'^l A = ( ^ 

We state a useful formula 

(25.7) ^(adj A) = (adj A) A = (det A)I. 

The (i, j) entry of A(adj A) is in fact the inner product of the i-th row of A and the i-th row of cof A, 
and thus (25.7) is simply the Laplace formula. 

We define the characteristic polynomial 

(25.8) pa{\) = det(yl - A/) 
of the matrix A. Let 

p^(A) = (-l)"A"+piA"-i + ---+Po. 
It naturally leads to a matrix polynomial 

An important property of the adjugate is 

(25.9) adj A = -((-1)"^'^-^ + piA""-^ + ■■■+ pil). 
Formally, adj A = ^^(t^. 
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